Absolute Measurement of Hadronic Branching Fractions of the D'^ Meson 
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The branching fractions of £)* meson decays serve to normahze many measurements of processes involving 
charm quarks. Using 298 pb"' of e^e" coUisions recorded at a center of mass energy of 4. 17 GeV, we determine 
absolute branching fractions for eight £)* decays with a double tag technique. In particular we determine the 
branching fraction S{D* — » K K*n^) = (5.50 ± 0.23 ± 0.16)%, where the uncertainties are statistical and 
systematic respectively. We also provide partial branching fractions for kinematic subsets of the K^K^n* decay 
mode. 

PACS numbers: 13.25.Ft 



Uncertainties in the decay probabilities (branching frac- 
tions) of the Dj meson to various detectable final states sig- 
nificantly impact the precision of a diverse array of measure- 
ments, including tests of the Standard Model prediction of the 
coupling of the Z" boson to charm quarks, measurements of 
B meson properties such as Z?'' mixing parameters, tests of 
light quark SU(3) symmetry in the D system, and tests of lat- 
tice gauge theory in leptonic decays. Any rate measure- 
ment where a meson is an intermediate step in a decay 



chain demands that the relevant normalizing branching frac- 
tions be known precisely to reduce systematic uncertainties. 
Most Z)^ branching fractions are presently obtained by com- 
bining measurements of ratios with a single absolute branch- 
ing fraction of one decay mode, thus introducing strong cor- 
relations and compounding uncertainties. In this Letter we 
present the first simultaneous high-statistics determination of 
multiple Dj absolute branching fractions, using a technique 
with significantly different systematic uncertainties from pre- 
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vious branching fraction measurements, which results in pre- 
cision better than current world averages by a factor of two. 
The eight decays considered in this analysis are — > /T" K^, 
K K+n^, K-K^n+7T^\ K'-^ K-n+n^, n+n+n-, n+ri, n^r]', and 
K^n^ji^ . Except where noted, mention of a decay implies 
the charge conjugate process as well. 

The most precise measurements of absolute branching 
fractions are currently obtained using partial reconstruction 
techniques to obtain the total number of D, mesons produced, 
either from B D^*}p[*^^^ decays Jll |2 or from e+e^ ^ 
D*='Djj(2536) events [3].'References and H quote results 
for the resonant decay (pn'^, while Ref. [30 measures 

Here we employ a technique extensively used by CLEO-c, 
pioneered by the MARK III collaboration for measuring 
and branching fractions P, ^ and limiting Dg branch- 
ing fractions |6], which exploits a feature of near-threshold 
production of charmed mesons. Below the DsDK thresh- 
old of 4.33 GeV, production of a meson in a charm- and 
strangeness-conserving process requires the existence of a Dj 
meson elsewhere in the event (possibly with additional pho- 
tons or pions). Events where at least one candidate is re- 
constructed ("single tag" or ST events) thus provide a sample 
with a known number of Ds events. Absolute branching frac- 
tions can then be obtained by finding events with two recon- 
structed Ds candidates ("double tag" or DT events). In this 
analysis, yields for charge-conjugate ST modes are consid- 
ered separately, but charge-conjugate branching fractions are 
assumed to be equal, ignoring the possibility of direct CP vi- 
olation. There are 16 ST yields, corresponding to two charges 
for each considered Ds decay, and 64 DT yields, one for each 
pairing of a D^ and a decay. 

This analysis uses (298 + 3) pb"' of data taken at a center 
of mass energy of 4.17 GeV. At this energy the dominant Ds 
production mechanism is the process e^e" D*s^D^s with 
a cross-section of ~ 1 nb f/J; the D* then decays to either 
yDs or n'^Ds in a ~ 16 : 1 ratio The very small rate of 
e^e^ — > D^Ds is not used for this analysis. The transition 
photon or n^^ is not reconstructed. 

To illustrate the method, consider two ST modes, D^ — > / 
and Z)j J, and one DT mode Z)^ — > /, D7 /. For a given 
number of produced D, pairs Nd^.d,, branching fractions S, 
and Sj, and efficiencies for the ST (e, and ej) and DT (e,|) 
events, we expect to observe the yields 



yj = No-DSjej, 



where y, and yj are the ST yields and y,,- is the DT yield. Using 
e,, fj, and eq from Monte Carlo simulations, we can solve for 



the branching fractions and Nd',d,' 

yj 

yi eij 

yiyj e, 



Nn 



yi] ^i^j 



In practice, to maximize the statistical power of the analysis, 
the parameters Nd',d, and S; are simultaneously optimized us- 
ing a maximum likelihood fit to the observed yields, where the 
ST yields use Gaussian likelihood functions and the consider- 
ably smaller DT yields use Poisson likelihood functions. The 
statistical properties of the fit were checked with pseudoex- 
periments. 

The CLEO-c detector is a modification of the CLEO III de- 
tector 1 iL II2II . The silicon strip vertex detector was 



replaced by a six layer vertex drift chamber. The charged par- 
ticle tracking system, consisting of the vertex chamber and 
a 47-layer central drift chamber, operates in an axial 1 T 
magnetic field, and provides a momentum resolution o-plp ~ 
0.6% 2A. p — \ GeV/c for tracks traversing every layer. Pho- 
tons are detected in an electromagnetic calorimeter consisting 
of 7784 CsI(Tl) crystals, which achieves an energy resolution 
o-eIE ~ 5% at 100 MeV. Two particle identification (PID) 
systems are used to distinguish charged kaons and pions: the 
central drift chamber, which provides specific ionization mea- 
surements (dE/dx), and, surrounding this chamber, a cylindri- 
cal Ring Imaging Cherenkov (RICH) detector. The combined 
PID system achieves n"^ and K"^ identification efficiency ex- 
ceeding 85% with fake rates less than 5% over the kinematic 
range of interest li 1 3i1 . The detector response is modeled with a 
detailed GEANX-based Hill Monte Carlo simulation, with initial 
particle trajectories generated by evtgen jlSR and final state 
radiation produced by photos 116!]. The initial state radiation 
spectrum is modeled using cross-sections for D*Ds produc- 
tion at lower energies determined during a CLEO-c scan of 
this region 

Charged tracks are required to be well-reconstructed and, 
except for K'^^ daughters, to be consistent with originating 
at the interaction point. The initial selection requires track 
momentum > 50 MeV/c. Candidate K'^ and n"^ tracks are 
chosen using dE/dx and RICH information, using the same 
criteria as the CLEO-c D^/D^ absolute branching fraction 
analysis 11311 . Charged kaons must have momentum above 
125 MeV/c. Neutral kaon candidates are reconstructed in the 
— > TT^n^ decay. The two pions have no PID require- 
ments, and a vertex fit is done to allow for the K^^ flight dis- 
tance. The pion pair is required to have an invariant mass 
within 6.3 MeV/c^ of the nominal K'^^ mass. We form tt" and 
77 candidates using pairs of isolated electromagnetic showers, 
keeping combinations within three standard deviations of the 
nominal masses; for further use a kinematic fit constrains the 
candidates to the nominal mass. Candidate 77' mesons are re- 
constructed by combining rj candidates with n^n^ pairs; the 
pions are subject to the standard pion PID requirements, and 
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the reconstructed 77' mass must be within lOMeV/c^ of the 
nominal value. 

We use several samples of simulated events to obtain effi- 
ciencies, study background shapes, and cross-check the analy- 
sis. A "generic" decay models a physical D, decaying into any 
of its final states; the branching fractions and intermediate res- 
onant components used for various final states are motivated 
by Particle Data Group (PDG) averages |8|]. A "signal" decay 
is one in which the simulated D, always decays to a final state 
of interest, with the same ratio of resonant components as in 
generic decays. We obtain efficiencies from samples with ei- 
ther one signal and one generic decay (ST modes) or two sig- 
nal decays (DT modes). Backgrounds are investigated using a 
combined sample of generic D*\ D^, and decays with ap- 
propriate production mechanisms and rates at 4.17 GeV, and 
samples of e^e^ — » t^t", y\jj{2S), and light quarks. 

We identify candidates using their momenta and invari- 
ant masses. A candidate may either be the daughter of a D* 
(an "indirect" Ds) or be produced in the inital e^e^ D*Ds 
process (a "direct" Ds). Direct Ds candidates have fixed mo- 
mentum in the center of mass frame because they are pro- 
duced in a two-body process; indirect candidates have a 
momentum distribution smeared around this value due to the 
extra boost of the D* —> (y, 7i^)Ds decay. We define the recoil 
mass variable Mrec through 
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where (Eq^Po) is the e^e" center of mass four-vector, po^ 
is the measured momentum, and Ma, is the nominal D, 
mass. For direct candidates, M^c peaks at the D* mass 
of 2.112 GeV/c^; for indirect candidates, Mjec spreads 
roughly ±60MeV/c^ around this peak. For DT and most 
ST candidates, we require Mjec > 2.051 GeV/c^; this ac- 
cepts all kinematically allowed events. For three ST modes 
(K^K^n^rf', n'^n^n^, K^n^n^) tighter mode-dependent se- 
lections ofMiec > (2.099, 2.101, 2.099) GeV/c^, which are 
roughly 80% efficient for signal, are applied to improve the 
signal to background ratio. The M,ec requirement eliminates 
contributions from e^e^ D^D^ events as those occur in a 
narrow peak at M,ec = Mo^ ■ 

The Ds candidates are subject to mode-dependent vetoes 
to reduce structure in the background invariant mass spec- 
trum, mostly arising from copiously produced D*D* events. 
In all modes except K'^ and K^K^n^, all neutral and 
charged pions, including K^^ daughters, must have momen- 
tum above lOOMeV/c to eliminate the soft pions from D* 
decays. Reflections are reduced by vetoing candidates where 
certain daughter combinations are consistent with the or 
masses (for example, the K K^ pair in a K^K^n^ candi- 
date must not be consistent with a Cabibbo-suppressed D" de- 
cay). To remove contamination from K^^ decays in the modes 
TT^n'^TT^ and K^tt'^tt^, no n'^n^ combination may have a mass 
between 475 and 520 MeV/c^. 

For ST yield extraction, every event is allowed to contribute 
a maximum of one Ds candidate per mode and charge. If there 
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FIG. 1 : (Color online) Invariant masses of the £)+ candidates in data 
in ST modes. Charge-conjugate modes are combined. The fits for 
yields are shown as the dashed red lines (background component) 
and solid blue lines (signal plus background). The total ST yield is 
(30.9 ± 0.3) X 10' events. 



are multiple candidates, the one with Mrec closest to M^. is 
chosen. An unbinned maximum likelihood fit is then per- 
formed on the invariant mass spectrum of the candidates in 
each mode. The expected signal distribution is obtained from 
Monte Carlo simulations; in fits to data the Ds mass is allowed 
to float. Backgrounds are modeled with linear functions for all 
modes except K^K'^7T^7^'^' and ;7r^7r^7r", where quadratic func- 
tions are used. The same background shape is used for both 
charges in a given mode. The reconstructed candidate masses 
M(Ds) and ST yield fits are shown in Fig.[Tl Efficiencies for 
ST modes range from 5.3% to 51%. 

Double tag yields are extracted by defining a signal re- 
gion in the two-dimensional plane of the two Ds candidate 
masses, M(D^) vs. M{D^). Every event is allowed to con- 
tribute at most one DT candidate per possible final state; 
amongst multiple candidates, the combination with average 
mass M - (M(D'^) + M(D;))/2 closest to Mo, is chosen. The 
distribution of MiD') versus M(D'^) for all DT candidates, 
along with the signal and sideband regions, is shown in Fig.|2] 
Combinatoric background tends to have structure in M but be 
flat in the mass difference AM - M(D^) - M(D^y, in par- 
ticular simulations verify that the multiple candidate selection 
does not cause backgrounds to peak in AM. Both signal and 
sideband regions require |M - M^ | < 12 MeV/c^. The sig- 
nal region is |AM| < 30 MeV/c^, while the sideband region is 
50 < |AM| < 140 MeV/c^. Efficiencies for DT modes range 
from 0.3% to 38%. 

The Ds decay final states under consideration can often be 
reached through multiple intermediate resonances. For exam- 
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FIG. 2: (Color online) Masses of the Z)^ and candidates for all 
64 DT modes in data. The rectangles show the signal region (center) 
and two sideband regions (diagonally offset). There are 1089 events 
in the signal region and 339 events in the combined sideband regions. 



pie, in our Monte Carlo the final state K K'^tt^ is an incoher- 
ent mixture of K*°K^ (43%), 07r+ (38%), Zo(1430)°/:+ (8%), 
nonresonant production (7%) and /()(980)7r^ (4%). The recon- 
struction efficiency can depend significantly on which reso- 
nances are produced. Knowledge of the relative contributions 
of these intermediate states is incomplete. We compare in- 
variant mass distributions of pairs of Ds daughters in data and 
Monte Carlo, and use the resulting information on resonant 
structures to reweight the assumed intermediate state compo- 
nents. The resulting excursions in the efficiency are taken as 
systematic uncertainties. Where there is a significant com- 
ponent that cannot be explicitly assigned to any intermediate 
state, we find the worst-case variations between the dominant 
components. As an illustration, for K K^n^ we find that (pn'^ 
and K have very similar (and lowest) efficiencies, while 
the nonresonant component is 7% higher and the others lie be- 
tween these extremes. By selecting on the K K^ and K^n^ in- 
variant masses we ascribe 90% of reconstructed events to 07r^ 

or K K^; varying the assumed efficiency for the remaining 
events within the limits above changes the infeiTed average 
efficiency, leading to a systematic uncertainty of 1.5%. The 
uncertainties assigned vary from zero for the two-body final 

states to 6% for K^K^n^n'^ (where there is a large efficiency 

— *o 

difference between (pp and K K ). We also include uncer- 
tainties in the PDG 2007 fit values for Jj) (0.7%) and 
'B(r\' — » n^n^rj) (3.1 %), and correct for the difference between 
the PDG fit for S(/r^ — > tt^tt ) and the value used in geant. 

Systematic uncertainties for the simulation of track, K^^ , 
tt", and rj reconstruction and FID efficiencies are determined 



using partial versus full reconstruction of events in CLEO-c's 
ifr{2S) and (^(3770) datasets; the methods are shared with the 
D^/D^ branching fraction analysis [13]. Tracking efficien- 
cies are verified using i^(3770) — » DD events for n"^ and 
K'^, and using t//(2S) —> n^TT^J/ifi for tt'^. Good agreement 
is found, and an uncertainty of 0.3% per track is used, cor- 
related among all tracks, with an additional uncertainty of 
0.6% per kaon added in quadrature. Systematic effects in 
the PID efficiency are studied using i/'(3770) DD events; 
in general data has slightly lower efficiency than the simu- 
lations and corrections are applied. Because the corrections 
are momentum-dependent this is also affected by the uncer- 
tainty on the intermediate resonant states. The coiTections 
appHed range from (-0.2 + 0.2)% for jt+tj to (-3.7 + 1.4)% 
for K'K'^n'^TT^. Neutral kaon efficiencies are verified using 
DD events and the D^ K'^ mode; a systematic uncer- 
tainty of 1.9% per candidate is used. The 7r° efficiency is 
checked with i/'(25) — » ji^jtPjIiI/ decays, and the t] efficiency 
with \j/{2S ) — > rjJIil/ events. In both cases there are discrepan- 
cies between data and the simulation, and relative coiTections 
of (-3.9 + 2.0)% per ifi and (-5.7 + 4.0)% per 77 are appHed. 

The nominal signal lineshapes used in the ST yield fits are 
derived from the simulation, and the backgrounds are either 
linear or quadratic. We determine systematic uncertainties in 
the yields by relaxing each assumption separately: the mass 
resolution is allowed to vary by an overall scale factor, and 
the background is parameterized by a second-order polyno- 
mial if the nominal fit uses a linear one, or vice versa. The 
size of the resulting excursions vary from 0.2% (K^K^n^) to 
8.6% (K-K^n^n^'') for background shape and 0.1% (K^^ K"-) to 
10.3% {n^rf} for width. 

The efficiency for a reconstructed DT event to lie in the 
signal region depends on the mass resolutions for both candi- 
dates. EiTors in modeling the resolution will thus cause errors 
in the DT efficiency which are coiTelated with the ST signal 
lineshape uncertainties. To estimate this effect we use the best 
fit results from the ST width check to determine the changes 
expected in the DT efficiency. The difference due to each de- 
cay mode is taken as a systematic uncertainty competely cor- 
related with the corresponding ST uncertainty. The range of 
these effects is 0-8%. 

In addition, we consider mode-dependent systematic uncer- 
tainties arising from our modeling of average D*fis event mul- 
tiplicity and detector noise (0-3%), the final state radiation 
spectrum generated by photos (0.2-1 .2%), and our simulation 
of initial state radiation (0-0.8%). 

Peaking backgrounds in ST events are found to be negli- 
gible compared to the size of the background shape uncer- 
tainties. Very small crossfeeds (of order 0.5% or less) are 
expected between various DT modes and are included in the 
fit; peaking DT backgrounds from other sources mostly arise 
from D*D* reflections and are again found to be negligible. 

Systematic uncertainties are propagated to the final results 
by altering fit inputs (efficiencies and yields) with appropri- 
ate correlations and noting the variations in the results. The 
analysis was validated on a simulated generic sample of open 
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TABLE I: Branching fraction results from this analysis, world average branching fractions from the PDG 2007 fit 1 8], ratios of branching frac- 
tions to B(D* K^K^n*), and charge asymmetries ^cp- Uncertainties on CLEO-c measurements are statistical and systematic, respectively. 



Mode 



This Result S (%) PDG 2007 fit S (%) S/S(K-K^7i^) 



1.49 ±0.07 ±0.05 

5.50 ±0.23 ±0.16 
5.65 ± 0.29 ± 0.40 
1.64 ±0.10 ±0.07 
1.11 ±0.07 ±0.04 
1.58 ±0.11 ±0.18 
3.77 ± 0.25 ± 0.30 
0.69 ± 0.05 ± 0.03 



2.2: 
5.3 : 



0.4 
0.8 



2.7 ±0.7 
1.24 ±0.20 
2.16 ±0.30 

4.8 ±0.6 
0.67 ±0.13 



0.270 ± 0.009 ± 0.008 
1 

1.03 ±0.05 ±0.08 
0.298 ±0.014 ±0.011 
0.202 ±0.011 ±0.009 
0.288 ±0.018 ±0.033 

0.69 ± 0.04 ± 0.06 



+4.9 ±2.1 ±0.9 
+0.3 ± 1.1 ±0.8 
-5.9 ± 4.2 ± 1.2 
-0.7 ± 3.6 ± 1.1 
+2.0 ± 4.6 ± 0.7 
-8.2 ±5.2 ±0.8 
-5.5 ± 3.7 ± 1.2 



0.125 ±0.009 ±0.005 +11.2 ± 7.0 ± 0.9 



TABLE II: Partial branching fractions Sam for K^K*n* events with 
K K* mass within AM MeV/c^ of the mass. Uncertainties are 
statistical and systematic, respectively. 



Value This Result S (%) 
'Ws 1.69 ±0.08 ±0.06 
S,o 1.99 ±0.10 ±0.05 
Si5 2.14 ±0.10 ±0.05 
S20 2.24 ±0.11 ±0.06 



charm production with thirty times the statistics of the data, 
and successfully reproduced the input branching fractions. 

We have separate yields and efficiencies for and 
events, so it is possible to compute asymmetries 



yi/ei +yilei' 



which are sensitive to direct CP violation in decays (ex- 
pected to be very small in the Standard Model). Most system- 
atic uncertainties cancel in this ratio; the ones that remain are 
due to charge dependence in tracking and PID, and the depen- 
dence of the ST yields on the signal lineshape and background 
parametrization. 

The obtained branching fractions, branching ratios, and CP 
asymmetries are shown in Table U The values we obtain 
are consistent with the world averages [8] and significantly 
more precise than any previous absolute measurements of 
branching fractions. This is also the first result where all eight 
modes are measured simultaneously; the PDG fit combines 
many disparate branching ratio results. No significant CP 
asymmetries are observed. We additionally obtain the num- 
ber of D*Ds events Nd-d, = (2.93 + 0.14 + 0.06) x 10^ which 
gives o-d;o,(4.17 GeV) = (0.983 + 0.046 + 0.021 + 0.010) nb; 
in order, the uncertainties are statistical, systematic due to this 
measurement, and for the cross-section, systematic due to lu- 
minosity measurement [13]. The cross-section is consistent 
with earlier CLEO-c results obtained via a scan of this energy 
region [7]. 

A quantity conventionally termed S(Z)^ <pJ^^) has often 
been used as a reference branching fraction for decays; 
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FIG. 3: (Color online) Yields of £»; -> K-Ji* single tag events 
versus K K* invariant mass; no efficiency corrections have been 
applied. The ST fit procedure for the full K K^n^ sample is ap- 
plied here to the subsample of each bin of M(K^K^) and the result- 
ing yields plotted, hence backgrounds have been subtracted and the 
yields shown are signal. A ip peak is visible above an additional broad 
signal component. The lines show the mass window boundaries for 
the partial branching fractions in Table llll 



operationally it is measured by making kinematic selections 
on the kaon pair in Dl. — » K K^n* events and assuming a 
pure K K^ signal. However, the Dalitz plot for this 
mode shows the presence of a significant broad scalar com- 
ponent under the peak, whose contribution to the observed 
yield varies from less than 5% to over 10% depending on the 
4> candidate selection criteria. Figure [3] shows the mass spec- 
trum of Dj K^K^TT^ events in this mass region; when fit 
by a single Gaussian, the M{K^K^) resolution is 1.1 MeV/c. 
The scalar component will additionally interfere with the 
contribution, altering the observed rate of events in the </> peak 
from the Dl. — > (f)7T^ fit fraction which would be measured 
in an amplitude analysis. These variations are comparable to 
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or exceed the systematic uncertainties in our measurements. 
For this reason, we do not quote a branching fraction for the 
resonant mode (pn^; this quantity can only be unam- 

biguously measured with an amplitude analysis, which is of 
limited utility as a reference branching fraction. We instead 
provide partial branching fractions Sam, which are defined as 
the branching fraction for K K^n^ events where the K 
pair satisfies \M(K-K+) - 1019.5 MeV/c^l < AM (MeV/c^); 
the values obtained are listed in Table HIl The systematic un- 
certainties quoted for Sam include uncertainties due to reso- 
lution. We emphasize that these are not measurements of the 
quantity S(D^ — > (pn^ —> K'K^k^), but are intended as ref- 
erences to normalize other Dj branching fractions when most 
of the K K^n^ phase space must be excluded for background 
reasons. 

In summary, we have presented the first application of a tag- 
ging technique at a center of mass energy of 4. 17 GeV to mea- 
sure eight hadronic branching fractions with precisions ex- 
ceeding world averages. For the key mode Dj — » K K^n^, 
the statistical and systematic uncertainties are comparable. As 
the experimental meaning of S(Z)^ — > (pn^) is ill-defined at 
this level of precision without a full amplitude analysis, we do 
not report it. We instead provide partial branching fractions 
for windows centered on the <p mass which do not assume a 
specific resonant composition of the decay. 
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